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O . Abstract 

! We consider a 1^ degrees of freedom Hamiltonian dynamical system, which models 

the chaotic dynamics of charged test-particles in a turbulent electric field, across the 
■ confining magnetic field in controlled thermonuclear fusion devices. The external electric 

field E = — W is modeled by a phenomenological potential V and the magnetic field B is 
considered uniform. It is shown that, by introducing a small additive control term to the 
external electric field, it is possible to create a transport barrier for this dynamical system. 
The robustness of this control method is also investigated. This theoretical study indicates 



Q^l that alternative transport barriers can be triggered without requiring a control action on the 



o 



device scale as in present Internal Transport Barriers (ITB). 



?3 : Contents 

1 Introduction 

QQ ■ 2 Physical motivations and the E x B model 2 

■ 2.1 Physical motivations 2 

^ : 2.2 The Ex B model 



Q 3 Localized control theory of hamiltonian systems 5 

. 3.1 The control term 5 

3.2 Properties of the control term |9 



^ ■ 4 Numerical investigations for the control term 

^ ! 4.1 Phase portrait for the exact control term 

4.2 Robustness of the barrier 

4.3 Energetical cost 

5 Discussion and Conclusion 

5.1 Main results 

5.2 Discussion, open questions M 

5.2.1 Comparison with the global control method 15| 

5.2.2 Effectiveness and properties of the control procedure Ts] 



*Centre de Physique Thcoriquc (CPT) - CNRS, Luminy - Case 907, 13288 Marseille cedex 09, France. 
Unite Mixte de Recherche (UMR 6207) du CNRS, et des universites Aix-Marseille I, Aix-Marseille II et du Sud 
Toulon- Var. Laboratoire affilic a la FRUMAM (FR 2291). Laboratoire de Recherche Conventionne du CEA 
(DSM-06-35). 

■^Institut dc Recherche sur la Fusion Magnctiquc (IRFM/SCCP), CEA Cadarachc, 13108 St-Paul-lcz-Durance, 
France. 

*M2P2, UMR 6181, Tcchnopole dc Chateau Gombcrt, 13451 Marseille Cedex 20, France. 



1 



1 Introduction 



It has long been recognised that the confinement properties of high performance plasmas with 
magnetic confinement are governed by electromagnetic turbulence that develops in microscales [1] . 
In that framework various scenarios are explored to lower the turbulent transport and therefore 
improve the overall performance of a given device. The aim of such a research activity is two-fold. 
First, an improvement with respect to the basic turbulent scenario, the so-called L-mode (L for 
low) allows one to reduce the reactor size to achieve a given fusion power and to improve the 
economical attractiveness of fusion energy production. This line of thought has been privileged 
for ITER that considers the H-mode (H for high) to achieve an energy amplification factor of 10 
in its reference scenario [2]. The H-mode scenario is based on a local reduction of the turbulent 
transport in a narrow regime in the vicinity of the outermoster confinement surface [3]. 
Second, in the so-called advanced tokamak scenarios. Internal Transport Barriers are considered 
[2] . These barriers are characterised by a local reduction of turbulent transport with two important 
consequences, first an improvement of the core fusion performance, second the generation of 
bootstrap current that provides a means to generate the required plasma current in regime with 
strong gradients [1]. The research on ITB then appears to be important in the quest of steady 
state operation of fusion reactors, an issue that also has important consequences for the operation 
of fusion reactors. 

While the H-mode appears as a spontaneous bifurcation of turbulent transport properties in the 
edge plasma [5], the ITB scenarios are more difficult to generate in a controlled fashion [3]. Indeed, 
they appear to be based on macroscopic modifications of the confinement properties that are both 
difficult to drive and difficult to control in order to optimise the performance. 
In this paper, we propose an alternative approach to transport barriers based on a macroscopic 
control of the E x B turbulence. Our theoretical study is based on a localized hamiltonian control 
method that is well suited for E x B transport. In a previous approach [6], a more global scheme 
was proposed with a reduction of turbulent transport at each point of the phase space. In the 
present work, we derive an exact expression to govern a local control at a chosen position in 
phase space. In principle, such an approach allows one to generate the required transport barriers 
in the regions of interest without enforcing large modification of the confinement properties to 
achieve an ITB formation [S] . Although the application of such a precise control scheme remains 
to be assessed, our approach shows that local control transport barriers can be generated without 
requiring macroscopic changes of the plasma properties to trigger such barriers. The scope of 
the present work is the theoretical demonstration of the control scheme and consequently the 
possibility of generating transport barriers based on more specific control schemes than envisaged 
in present advanced scenarios. 

In Section [21 we give the general description of our model and the physical motivations for our 
investigation. In Section [3l we explain the general method of localized control for Hamiltonian 
systems and we estimate the size of the control term. Section H] is devoted to the numerical 
investigations of the control term, and we discuss its robustness and its energy cost. The last 
section [5] is devoted to conclusions and discussion. 

2 Physical motivations and the E x B model 
2.1 Physical motivations 

Fusion plasma are sophisticated systems that combine the intrinsic complexity of neutral fiuid 
turbulence and the self-consistent response of charged species, both electrons and ions, to magnetic 
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fields. Regarding magnetic confinement in a tokamak, a large external magnetic field and a 
first order induced magnetic field are organised to generate the so-called magnetic equilibrium 
of nested toroidal magnetic surfaces [7j. On the latter, the plasma can be sustained close 
to a local thermodynamical equilibrium. In order to analyse turbulent transport we consider 
plasma perturbations of this class of solutions with no evolution of the magnetic equilibrium, 
thus excluding MHD instabilities. Such perturbations self-consistently generate electromagnetic 
perturbations that feedback on the plasma evolution. Following present experimental evidence, 
we shall assume here that magnetic fluctuations have a negligible impact on turbulent transport 
[8]. We will thus concentrate on electrostatic perturbations that correspond to the vanishing (5 
limit, where (3 = p/(i?^/2/io) is the ratio of the plasma pressure p to the magnetic pressure. The 
appropriate framework for this turbulence is the Vlasov equation in the gyrokinetic approximation 
associated to the Maxwell-Gauss equation that relates the electric field to the charge density. 
When considering the Ion Temperature Gradient instability [9] that appears to dominate the ion 
heat transport, one can further assume the electron response to be adiabatic so that the plasma 
response is governed by the gyrokinetic Vlasov equation for the ion species. 
Let us now consider the linear response of such a distribution function /, to a given electrostatic 
perturbation, typically of the form Tg ^-ti^t+ikr ^ (where / and are Fourier amplitudes of 
distribution function and electric potential). To leading orders one then finds that the plasma 
response exhibits a resonance: 



Here feq is the reference distribution function, locally Maxwellian with respect to f|| and uj* is 
the diamagnetic frequency that contains the density and temperature gradient that drive the ITG 
instability [9]. Tg is the electronic temperature. This simplified plasma response to the electrostatic 
perturbation allows one to illustrate the turbulent control that is considered to trigger off transport 
barriers in present tokamak experiments. 

Let us examine the resonance u — k\\ v\\ = where k\\ = {n — m/q)/R with R being the major 
radius, q the safety factor that characterises the specific magnetic equilibrium and m and n the 
wave numbers of the perturbation that yield the wave vectors of the perturbation in the two 
periodic directions of the tokamak equilibrium. When the turbulent frequency uj is small with 
respect to Vth/{qR), (where Vth = \/kBT jva is the thermal velocity), the resonance occurs for 
vanishing values of k\\^ and as a consequence at given radial location due to the radial dependence 
of the safety factor. The resonant effect is sketched on figure 1. 



Figure 1: Resonances for q = — and q = -^^d— for two different widths, narrow resonances 
empedding large scale turbulent transport and broad resonances favouring strong turbulent 
transport. 

In a quasilinear approach, the response to the perturbations will lead to large scale turbulent 
transport when the width of the resonance 6m is comparable to the distance between the resonances 






^m,m+i leading to an overlap criterion that is comparable to the well known Chirikov criterion 
for chaotic transport am. = (^m + ^m+i) / ^m,m+i with a > 1 leading to turbulent transport across 
the magnetic surfaces and a < 1 localising the turbulent transport to narrow radial regions in the 
vicinity of the resonant magnetic surfaces. 

The present control schemes are two- fold. First, one can consider a large scale radial electric field 
that governs a Doppler shift of the mode frequency u. As such the Doppler shift u — uoe has no 
effect. However a shear of the Doppler frequency uje, ^jJe = uJE + ^foj'^ will induce a shearing effect 
of the turbulent eddies and thus control the radial extent of the mode so that one can locally 
achieve a < 1 in order to drive a transport barrier. 

Second, one can modify the magnetic equilibrium so that the distance between the resonant 
surfaces is strongly increased in particular in a magnetic configuration with weak magnetic shear 
{dq/dr ^ 0) so that A^.m+i is strongly increased, Am,m+i ^ (^m, also leading to a < 1. 
Both control schemes for the generation of ITBs can be interpreted using the situation sketched 
on figure 1. The initial situation with large scale radial transport across the magnetic surfaces (so 
called L-mode) is indicated by the dashed lines and is governed by significant overlap between the 
resonances. The ITB control scheme aims at either reducing the width of the islands or increasing 
the distance between the resonances yielding a situation sketeched by the plain line in figure 1 
where the overlap is too small and a region with vanishing turbulent transport, the ITB, develops 
between the resonances. 

Experimental strategies in advanced scenarios comprising Internal Transport Barriers are based on 
means to enforce these two control schemes. In both cases they aim at modifying macroscopically 
the discharge conditions to fulfill locally the a < 1 criterion. It thus appears interesting to devise 
a control scheme based on a less intrusive action that would allow one to modify the chaotic 
transport locally by the choice of an appropriate electrostatic perturbation hence leading to a 
local transport barrier. 



For fusion plasmas, the magnetic field B is slowly variable with respect to the inverse of the Larmor 
radius pi i.e: pi|Vlni?| <C 1. This fact allows the separation of the motion of a charged test 
particle into a slow motion (parallel to the lines of the magnetic field) and a fast motion (Larmor 
rotation). This fast motion is named gyromotion, around some gyrocenter. In first approximation 
the averaging of the gyromotion over the gyroangle gives the approximate trajectory of the charged 
particle. This averaging is the guiding-center approximation. 

In this approximation, the equations of motion of a charged test particle in the presence of a 
strong uniform magnetic field B = Bz, (where z is the unit vector in the z direction) and of an 
external time-dependent electric field E = — Wi are: 



where Vi is the electric potential. The spatial coordinates X and Y play the role of canonically- 
conjugate variables and the electric potential Vi{X, Y, T) is the Hamiltonian for the problem. Now 
the problem is placed into a parallelepipedic box with dimensions L x i x {2tt/uj), where L and 



2.2 The E X B model 
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^ are some characteristic lengths and is a characteristic frequency of our problem, X is locally 
a radial coordinate and F is a poloidal coordinate. A phenomenological model [TU] is chosen for 
the potential: 

vi(A-.r,r)^E|^ (3) 

n,m=l ^ ' 

where Vq is some amplitude of the potential, 

uj is constant, for simplifying the numerical simulations and 0n,m are some random phases 

(uniformly distributed) . 

We introduce the dimensionless variables 

(x, y, t) = {2-kX/L, ujT) (4) 

So the equations of motion ([2]) in these variables are: 



— ( ^ \ = ( -dy^i^^y^'t) 



(5) 



where V = e{Vi/Vo) is a dimensionless electric potential given by 

cos {nx + my + (j)n,m - t) 



yr-^ (:ubyux ^ iiiy ^ qjn,m — 1^ 
V{x,y,t) = e > 6 

n,m=l ^ ' 

Here 

e = ATr\cVo/B)/{L£u) (7) 

is the small dimensionless parameter of our problem. We perturb the model potential ([6]) in order 
to build a transport barrier. The system modeled by Eqs.(l5]) is a l| degrees of freedom system 
with a chaotic dynamics [101 E]- The poloidal section of our modeled tokamak is a Poincare section 
for this problem and the stroboscopic period will be chosen to be 27r, in term of the dimensionless 
variable t. 

The particular choice ([3]) or ([6]) is not crucial and can be generalized. Generally, u can be chosen 
depending on n, m. This would make the numerical computations more involved. In the following 
section, V is chosen completely arbitrary. 



3 Localized control theory of hamiltonian systems 
3.1 The control term 

In this section we show how to construct a transport barrier for any electric potential V. The 
electric potential V{x, y, t) yields a non-autonomous Hamiltonian. We expand the two-dimensional 
phase space by including the canonically-conjugate variables {E,t), 

H = H{E,x,y,T) = E + V{x,y,T) (8) 



The Hamiltonian of our system thus becomes autonomous. Here r is a new variable whose 
dynamics is trivial: f = 1 i.e. r = tq + 1 and E is the variable canonically conjugate to r. The 
Poisson bracket in the expanded phase space for any W = W{E, x, y, r) is given by the expression: 

{W} = {d,W)dy - idyW)d, + idEW)dr - idrW)dE. (9) 

Hence {W} is a linear (differential) operator acting on functions of {E,x,y,T). We call Hq = E 
the unperturbed Hamiltonian and V{x,y,T) its perturbation. We now implement a perturbation 
theory for Hq. The bracket (Q for the Hamiltonian H is 

{H} = {d,V)dy - {dyV)d, + dr- {drV)dE (10) 

So the equations of motion in the expanded phase space are: 

y = {H}y = d,V{x,y,T) (11) 

X = {H}x = -dyVix,y,T) (12) 

E = {H}E = -d^V{x,y,T) (13) 

i = {H}t = 1 (14) 

We want to construct a small modification F of the potential V such that 

H = E + V{x, y, t) + F{x, y,r) = E + V{x, y, r) (15) 

has a barrier at some chosen position x = Xq. So the control term 

F = V{x,y,T)-V{x,y,T) (16) 

must be much smaller than the perturbation (e.g., quadratic in V). One of the possibilities is: 

V = V{x + dyf{y,T),y,T) (17) 

where 



/(y,^) = / V{xo,y,t)dt 
Jo 

Indeed we have the following theorem: 

Theorem 1 The Hamiltonian H has a trajectory x = Xq + dyf{y, r) acting as a barrier in phase 
space. 

Proof 

Let the Hamiltonian H = exp{{f})H be canonically related to H. (Indeed the exponential of any 
Poisson bracket is a canonical transformation.) We show that H has a simple barrier a.t x = Xq. 
We start with the computation of the bracket for the function /. Since / = f{y,T), the 
expression for this bracket contains only two terms, 

{/} ^ -fd. - fdE (18) 

where 

/' = dyf and / = d^f (19) 

which commute: 

[fd,JdE] = Q (20) 



Now let us compute the coordinate transformation generated by exp ({/}): 

exp({/}) = exp(-/'9^) exp{-fdE), 



(21) 



where we used ( 12U]) to separate the two exponentials. 

Using the fact that exp{bdx) is the translation operator of the variable x by the quantity b: 
[exp{bdx)W]{x) =W{x + b),we obtain 

H = e^f^H = e^f^E + e^f^V{x,y,T) 
E - f) + V{x-f',y,r) 
E - y(xo,y,r) + V{x + f'-f',y,T) 



E-V{xo,y,T) + V{x,y,T) 



(22) 



This Hamiltonian has a simple trajectory x = Xo,E = Eq, i.e. any initial data x = Xq, y = 
yo, E = Eq, t = Tq evolves under the flow of H into x = Xo,y = yt, E = Eq, r = Tq + 1 for some 
evolution yt that may be complicated, but not useful for our problem. Hamilton's equations for x 
and E are now 



X = {H}x = dy[V{xo,y,r) -V{x,y,T)] 
E = {H}E = dr[V{xQ,y,T)-V{x,y,T)] 



(23) 
(24) 



so that for x = xq, we find x = = E. Then the union of all points (x, y, E , t) at x = xq E = Eq: 

/ ^0 \ 



«o= U 



y 

Eq 

\ ^ I 



(25) 



is a 3-dimensional surface x M, (T = M/27rZ) preserved by the flow of H in the 4-dimensional 
phase space. If an initial condition starts on *Bo, its evolution under the flow exp(t{ff}) will 
remain on ^Bq. 

So we can say that *Bo act as a barrier for the Hamiltonian H: the initial conditions starting inside 
QSo can't evolve outside *Bo and vice-versa. ^ 

To obtain the expression for a barrier 03 for H we deform the barrier for B. via the transformation 
exp({/}). As 

E = e-^f^H (26) 
and exp({/}) is a canonical transformation, we have 



{H} = {e~^f^H} = e^^f^{H}e^f^ 
Now let us calculate the flow of H: 



AH} ^ t{e'yHH}eifi) ^ -{f}J{H} {f} 



Indeed: 



^t{e-if){H}e^n) ^ J2 

n=0 



t^{e-if}{H}e^f})^ 



(27) 

(28) 
(29) 
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For instance when n = 2: 



{/} 



{/} 



and so 



As we have seen before: 



ot{H} 



n=0 



,{/} 



and 



AH} 



( x\ 

y 

E 
\ ^ J 

y 

V ^ / 



Multiplying fl28l) on the right by e we obtain: 

/ ^0 \ 



y 

E-f 

V r ) 



yt 

Eo 
\r + t / 



and 



e-if}^t{H} 



y 

Eo 

V ^ / 

f Xo \ 

y 

Eo 

\ ^ ) 



Am 



( Xo + /'(y,r) \ 

y 

Eo + fiy,r) 



f Xo \ 



-{/} 



yt 

Eo 
\r + t / 

( xo + f'iyuT + t) \ 

yt 

Eo + f{yt, r + t) 

T + t 



So the flow exp{t{H}) preserves the set 



/ x„ + r{y,T) \ 



s= U 



y 

Eo + fiy,r) 

T 



J 



(30) 
(31) 



(32) 



(33) 



(34) 



(35) 



03 is a 3 dimensional invariant surface, topologically equivalent to x M into the 4 dimensional 
phase space. 03 separates the phase space into 2 parts, and is a barrier between its interior and 
its exterior. ^ is given by the deformation exp({/}) of the simple barrier ^o- 
The section of this barrier on the sub space {x,y,t) is topologically equivalent to a torus T^. 



This method of control has been successfully applied to a real machine: a traveling wave tube to 
reduce its chaos [TT] . 
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3.2 Properties of the control term 

In this Section, we estimate the size and the regularity of the control term (fT6|) . 



Theorem 2 For the phenomenological potential ^ the control term fTS\} verifies: 



\\F\\± ± < e^N^ — 



N ' N 47r 



(36) 




F = V{x + f')-V{x) = [ dsd,V{x + sf',y,T)f'{y,T) 



(37) 



and then use Cauchy's Theorem. 

4 Numerical investigations for the control term 

In this Section, we present the results of our numerical investigations for the control term F. 
The theoretical estimate presented in the previous section shows that its size is quadratic in the 
perturbation. Figure [2] shows the contour plot of V{x,y,t) and V{x,y,t) {V = V + F) ed, some 
fixed time t, for example ^ = f • One can see that the contours of both potentials are very similar. 
But the dynamics of the systems with V and V are very different. 

For all numerical simulations we choose the number of modes = 25 in ([6]). In all plots the 
abscissa is x and the ordinate is y. 




4.1 Phase portrait for the exact control term 

To explore the effectiveness of the barrier, we plot (in Fig. [3]) the phase portraits for the original 
system (without control term) and for the system with the exact control term F. We choose 
the same initial conditions. The time of integration is T = 2000, the number of trajectories: 



Figure 2: Uncontrolled and controlled potential for e = 0.6, t = f 



, xo = 2. 



Q 



Ntraj = 200 (number of initial conditions, all taken in the strip — 1 — vr < x < — vr; < y < 27i) 
and the parameter e = 0.9. We choose the barrier at position xq = 2. And to get a Poincare 
section, we plot the poloidal section when t G 27rZ. Then we compare the number of trajectories 
passing through the barrier during this time of integration for each system. We eliminate the 
points after the crossing. For the uncontrolled system 68% of the initial conditions cross the 
barrier at xq = 2 and for the controlled system only 1% of the trajectories escape from the zone 
of confinement. The theory announces the existence of an exact barrier for the controlled system: 
these escaped trajectories (1%) are due to numerical errors in the integration. 





Figure 3: Phase portraits without control term and with the exact control term, for e = 0.9, 

Xo = 2, Ntraj = 200. 

One can observe that the barrier for the controlled system is a straight line. In fact this barrier 
moves, its expression depends on time: 



X 



xo + f'iy,t) 



(3J 



But when t G 27rZ its oscillation around x = xq vanishes: f'{y,2k7r) = J^^^ dyV{xo,y,t)dt = 0. 
This is what we see on this phase portrait. In fact we create 2 barriers at position x = xq, and 
X = xo — 27T (and also at xq + 2mi) because of the periodicity of the problem. We note that the 
mixing increases inside the two barriers. The same phenomenon was also observed in the control 
of fluids [13], where the same method was applied. 



4.2 Robustness of the barrier 

In a real Tokamak, it is impossible to know an analytical expression for electric potential V. So we 
can't implement the exact expression for F. Hence we need to test the robustness of the barrier 
by truncating the Fourier decomposition (for instance in time) of the controlled potential. 



Fourier decomposition 

Theorem 3 The potential can be decomposed as V = Ylk&'^k, where 



N 



V, = e 



Jk{np) 



n,m=l ^ ' 



COS (?7 + A;e + (A; - \)t) 



(39) 



in 



with 



Vn,m{y) = nx + my + (f)n,m + neFc (40) 



= E 7 (41) 

n,m=l 

^'.(2') = E 7 f::"^c^' (42) 

n,m=l ^ ' 

Km,n,y = TLXq + my + (j)n,m (43) 

and J'k is the Bessel's function 

1 r 

Jk{np) = — I cos {ku — np sin u) du (44) 
Jo 



Proof We rewrite explicitly the expression f[T7jl for our phenomenological controlled potential 
Vix,y,ty. 



^ cos {n{x + f'{y,t)) + my + (pn.,' 



t 



V{x,y,t) =ey] ^ ^ (45) 

n,m=l ^ ' 

with 



m (cosK„,„,j^ - cos(i^„,^,y - t) ) 

/'(f. «) = ^ E ^ — pt;;;^^ («) 

n,m=l ^ ' 

With the definition (jH]) and (j42]) we have: 

/'(2/,t) = £(F,(2/) (1 - cost) - F,(2/) sint) (47) 

Let us introduce 

p = e{Fl^Fl)'l^ (48) 

and 6 by 

psinQ = -eF,{y) pcose = -eF,{y) (49) 

so that 

" eo.(,,-i + np.in(e + t)) 

^ (n^ + m2)3/2 ^ ^ 

n,m=l ^ ' 



Using Bessel's functions properties [T] 

cos(p sin 9) = J4(p) cos /c0 (51) 

sm{p sin O) = ^ Jk{p) sin A;e (52) 

we get 

cos [t] — t + np sin(9 + t)) = Jki^p) cos (^) (53) 

fcez 

where ^ = ?7 + A;0 + (A; — l)t, and we finally obtain (l39l) . The theorem is proved. ■ 
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During numerical simulations we truncate the controlled potential by keeping only its first 3 
temporal Fourier's harmonics: 



N 



Vtr=eY, 



Aq + Ai cos t + Bismt + A2 cos 2t + B2 sin 2t 



n,m=l 



m 



!)3/2 



(54) 



^0 = J(i{np) cos(?7 + 6) 

Ai = Joinp) cos 7] + J2{np) cos{t] + 26) 

Bi = Jo{np) sinr] - J2{np) ?,in{7] + 20) 

M = Jsinp) cos(?7 + 36) - Ji{np) cos(?7 - 6) 

B2 = -Jsinp) sin(?7 + 36) - Ji{np) sin(?7 - 6) 




0.05 



\ 
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Figure 4: Exact Control Term and Truncated Control Term with e 



0.6,t=f. 



Figure H] compares the two contour plots for the exact control term and the truncated control 
term (1541) . Figur^ compares the two phase portraits for the system without control term and for 
the system with the above truncated control term (15^ . The computation of Vtr on some grid has 
been performed in Matlab and the numerical integration of the trajectories was done in C. 




Figure 5: e = 0.3, T = 2000, Ntraj = 50. 
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Table 1: Squared ratios of the amplitudes of the control term and the uncontrolled electric potential 
CexXtr] ratios of electric energy of the control term and the uncontrolled electic potential rjf.^j'qtr] 
for the system with exact and truncated control term. 



e 




Ctr 




rjtr 


0.3 


0.1105 


0.1193 


0.6297 


0.1431 


0.4 


0.1466 


0.1583 


0.7145 


0.2393 


0.5 


0.1822 


0.1967 


0.8161 


0.3550 


0.6 


0.2345 


0.2137 


0.9336 


0.4883 


0.7 


0.2518 


0.2716 


1.0657 


0.6375 


0.8 


0.2858 


0.3038 


1.2119 


0.8014 


0.9 


0.3191 


0.3439 


1.3722 


0.9796 


1.5 


0.5052 


0.5427 


2.6247 


2.3037 



One can see a barrier for the system with the truncated control term. As for the system with the 
exact control term we create two barriers at positions x = xq and x = xq — 2tt and the phenomenon 
of increasing the mixing inside the barriers persist. 

4.3 Energetical cost 

As we have seen before, the introduction of the control term into the system can reduce and even 
stop the diffusion of the particles through the barrier. Now we estimate the energy cost of the 
control term F and the truncated control term Ftr = Vtr — V. 

Definition 1 The average of any function W = W{x, y, r) is defined by the formula: 

p2TT p2iT /'27r 

<\W\>= j dx dy dt\W{x,y,t)\ (55) 
Jo Jo Jo 

Now we calculate the ratio between the absolute value of the truncated control (electric potential) 
or the exact control and the uncontrolled electric potential: 

=< \F\^ > I < i^r > 

and 

Ctr =< \Ftr\^ > / < \V\^ > 

We also compute the ratio between the energy of the control electric field and the energy of the 
uncontrolled system in their exact and truncated version 

r]ex =< |VFp > / < |VVp > 

and 

ritr =< \VFtr\^ > I < |VVp > 
for different values of e. Results are shown in Table [H 

One can see that the truncated control term needs a smaller energy than the exact control term. 
In Table [21 we present the number of particles passing through the barrier in function of e, after 
the same integration time. 

Let AJ\f = Mwithout — -^tr be the difference between the number of particles passing through the 
barrier for the system without control and with the truncated control and A?7 = rjex — Vtr the 
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Table 2: Number of escaping particles without control term Afwuhout, and for the system with the 
exact control term Mexact and the truncated control term Mtr- 



e 


■^without 


^ exact 




0.4 


22% 


0% 


6% 


0.5 


26% 


0% 


18% 


0.9 


68% 


1% 


44% 


1.5 


72% 


1% 


54% 



Table 3: Difference ZW^ of the number of particles passing trough the barrier and difference of 
relative electric energy Ar/ for the controlled and uncontrolled system. 



e 




Ar] 


0.3 


8% 


0.49 


0.4 


16% 


0.47 


0.5 


8% 


0.46 


0.9 


24% 


0.39 


1.5 


18% 


0.32 



difference between the relative electric energy for the system with the exact control term and the 
system with the truncated control term. In Table [3] we present AJ\f and Ar; for differents values 
of e. 

For e below 0.2 the non controlled system is rather regular, there is no particles stream through 
the barrier, so we have no need to introduce the control electric field. For e between 0.3 and 0.9 
the truncated control field is quite efficient, it allows to drop the chaotic transport through the 
barrier by a factor 8% to 24% with respect to the uncontrolled system and it requires less energy 
than the exact control field. For e greater than 1 the truncated control field is less efficient than 
the exact one, because the dynamics of the system is very chaotic. For example when e = 1.5, 
there are 72% of the particles crossing the barrier for the uncontrolled system and 54% for the 
system with the truncated control field. At the same time the energetical cost of the truncated 
control field is above 70% of the exact one, which allows to stop the transport through the barrier. 
So for 5 > 1 we need to use the exact control field rather than the truncated one. 

5 Discussion and Conclusion 

In this article, we studied a possible improvement of the confinement properties of a magnetized 
fusion plasma. A transport barrier conception method is proposed as an alternative to presently 
achieved barriers such as the H-mode and the ITB scenarios. One can remark, that our method 
differs from an ITB construction. Indeed, in order to build-up a transport barrier, we do not 
require a hard modification of the system, such as a change in the q-profile. Rather, we propose 
a weak change of the system properties that allow a barrier to develop. However, our control 
scheme requires some knowledge and information relative to the turbulence at work, these having 
weak or no impact on the ITB scenarios. 

5.1 Main results 

First of all we have proved that the local control theory gives the possibility to construct a transport 
barrier at any chosen position x = Xq for any electric potential V{x,y,t). Indeed, the proof given 
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in section [3] does not depend on the model for the electric potential V. In Subsection 13.11 we 
give a rigorous estimate for the norm of the control term F, for some phenomenological model 
of the electric potential. The introduction of the exact control term into the system inhibits the 
particle transport through the barrier for any e while the implementation of a truncated control 
term reduces the particle transport significantly for e G (0.3, 1.0). 

5.2 Discussion, open questions 

5.2.1 Comparison with the global control method 

Let us now compare our approach with the global control method ||6j which aims at globally 
reducing the transport in every point of the phase space. Our approach aims at implementing 
a transport barrier. However, one also observes a global modification of the dynamics since the 
mixing properties seem to increase away from the barriers. 

Furthermore, in many cases, only the first few terms of the expansion of the global control term 
[6] can be computed explicitly. Here we have an explicit exact expression for the local control 
term. 

5.2.2 Effectiveness and properties of the control procedure 

In subsection \2.2\ we have introduced the dimensionless variables (jlj) and defined a dimensionless 
control parameter e = Aw'^ (cVq / B) / (Liu) . In the simplifying case where I = L = 211 /k is the 
characteristic length of our problem, we have e = ck^Vo/ {ujB). Let us consider a symmetric 
vortex, hence with characteristic scale 1/k. Let us now consider the motion of a particle governed 
by such a vortex. The order of magnitude of the drift velocity is therefore ve = kcVo/B and 
the associated characteristic time tett, tett = ^/{kvE), is the eddy turn over time. Let u be 
the characteristic evolution frquency of the turbulent eddies, here of the electric field, then the 
Kubo number K is K = 1/ujtett- This parameter is the dimensionless control parameter of this 
class of problems, and we remark that in our case K = e. It is also important to remark that 
the parameter K also characterises the diffusion properties of our system. Indeed, let 5 be a step 
size of our particle in a random walk process and let r be the associated characteristic time, the 
diffusion coefficient is then D = 6'^/t. Since one can relate the characteristic step and time by the 
velocity, 6 = VeT, on also finds: 

(vet)' Pc'Vi _J_ , 

We also introduce the reference diffusion coefficient D = /c~^a;, so that: 

D/D = K'^uT (57) 

They are two asymptotic regimes for our system. The first one, is the regime of weak turbulence, 
characterised by ujtett ^ 1 and therefore K <^ 1. In this regime, the electric potential evolution 
is fast, the particle trajectories only follow the eddy geometry on distances much smaller than 
the eddy size. The steps 6 are small and the characteristic time r of the random walk such that 
ur ~ 1. The particle diffusion (1571) is then such that: 

D/D ^ for UJTETT > 1 (58) 

The second asymptotic regime is the regime of strong turbulence, with uotett ^ 1 and K ^ 1. 
Particles then explore the eddies before decorrelation and the characteristic time of the random 
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step is typically r ~ tett and: 

D/D ^ K for ujtett < 1 (59) 

The first regime corresponds to the weak turbulence limit with weak Kubo number and particle 
diffusion and the second to strong turbulence and large Kubo number and particle diffusion. The 
control method developed in this article does not depend on K = e. There is always a possibility 
to construct an exact transport barrier. However for the numerical simulations, we have remarked, 
that for small e one can observe a stable barrier without escaping particles, and for e close or more 
than 1 there is some leaking of particles across the barrier. The barrier is more difficult to enforce. 
Also when considering the truncated control term, one finds that the control term is ineffective in 
the strong turbulence limit. 

Let us now consider the implementation of our method to turbulent plasmas where the turbulent 
electric field is consistent with the particle transport. The theoretical proof of an hamiltonian 
control concept is developped provided the system properties at work are completely known. For 
example the analytic expression for the electric potential. This is impossible in a real system, 
since the measurements take place on a finite spatio-temporal grid. This has motivated our 
investigation of the truncated control term by reducing the actually used information on the 
system. As pointed out previously, one finds that this approach is ineffective for strong turbulence. 
Another issue is the evolution of the turbulent electric field following the appearance of a transport 
barrier. This issue would deserve a specific analysis and very likely updating the control term on a 
trasnport characteristic time scale. An alternative to such a process would be to use a retroactive 
Hamiltonian approach (a classical field theory) [15] and to develop the control theory in that 
framework. 
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